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Abstract 

The properties of a vortex in a rotating superfluid Fermi gas are studied in the unitary limit. 
A phenomenological approach based on Ginzburg-Landau theory is developed for this purpose. 
The density profiles, including those of the normal fluid and superfluid, are obtained at various 
temperatures and rotation frequencies. The superfluid and normal fluid densities can be identified 
from the angular momentum density. The total free energy and angular momentum of the vortex 
are also obtained. 



*Electronic address: yinlan@pku.edu.cn 



1 



I. INTRODUCTION 



The field of BEC-BCS crossover has attracted great attention from both experimental 
and theoretical sides. In experiments, the system is studied near Feshbach resonances so that 
the scattering length can be tuned by an applied magnetic field. At low temperatures, when 
the scattering length is tuned from negative side to positive side, the system evolves from a 
BCS state to a molecular BEC state. The observation of molecular BEC was first reported 
in references 111 El- When the scattering length was fast tuned from the negative side to 

nn 

the positive side of the resonance, condensation of Fermion pairs was observed |3L|4|. Since 
then quite a few experiments have focused on variousproperties of the system, including 
collective excitations HQ, Sin «e exc itationS HB. and propel 
The observation of vortex lattice jl(| has provided a conclusive evidence of superfluidity 
across the resonance in this system. 
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The vortex structure in superfluid fermi gas have been studied theoretically jll . 

12, 3- However most of these works focus on the non-rotating case, which is different 



from the experimental situation. To consider the effect of rotation, we have developed a 
phenomenological approach, based on our previous work on the phase-slip phenomenon in 
the resonant Fermi superfluid 19 1. 

In this paper, we focus on the properties of a vortex in the unitary limit |20|, because it 
is a special strong-interaction region where the magnitude of the scattering length is much 
bigger than the interparticle distance. In the unitary limit, the interparticle distance replaces 
the scattering length as the most important length scale, and the physical properties become 
unitary. In the following, we first present our phenomenological approach and then study 
the vortex structure at various temperatures and rotation frequencies. 



II. FORMALISM 

In our previous work , we developed a phenomenological theory for the resonant Fermi 
superfluid based on Ginzburg-Landau theory. To consider the effect of rotation, we study 
the system in the rotating frame, in which the single-particle Hamiltonian is given by 

H — QL Z = ^-{p 2 ± + p 2 z ) + \rn{u\r 2 ± + u z z 2 ) -Qz-r ± xp ± , (1) 
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where = (p x ,p y , 0), = (x, y, 0), u± and u z are the transverse and longitudinal trapping 
frequencies, and Q is the rotation frequency. Here the trapping potential is assumed to have 
the rotation symmetry in the transverse plane. The single-particle Hamiltonian can be 
rewritten as 

H-nL z = J-(p - + l ~m[(ul - Q 2 y ± + u z z% (2) 

where p — eA/c is the canonical momentum, and eA/c = mQz x rj_. From Eq. (J2J), we 
conclude that the motion of a charge-neutral fermion in a rotating frame is analogous to an 
electron moving in a vector potential A and in a trapping potential reduced in the transverse 
direction. 

In the superfluid state, the free energy of a fermi gas is lower than that in the normal 
state, and the energy difference is defined as the condensation energy. The total free energy 
F is the sum of condensation energy Fq and the free energy of the normal state Fn, 

F = F C + F N . (3) 

The free energy of the normal state F^ is generally a complicated function of tempera- 
ture, particle density, trapping potential, the scattering length, and the rotation frequency. 
According to Ginzburg-Landau theory, the condensation energy can be expressed as a func- 
tional of the vector potential A and the order parameter <fi which is proportional to the 
energy gap of the fermion excitation, 

Fo = J rfM-|V(r)(V - + «(r)|0(r)| 2 + ^-\4>(r)\\ (4) 

where a and (3 are coefficients. Although Ginzburg-Landau theory was originally proposed 
for describing superconducting system near the transition, it can often be applied to most 
of the region below the transition temperature. 

When the system is not rotating, Q = 0, the sign change of the coefficient a signals the 
superfluid transition. Below the transition temperature T@, a < 0, and the system has lower 
free energy in the superfluid state than in the normal state. When the system is rotating, 
the transition temperature becomes a function of the rotation frequency Q. In the following, 
the symbol Tq is defined as the transition temperature at O = 0. 

In the superfluid state, the total free energy F is at a minimum point as a functional of 
the particle density n(r) and the superfluid order parameter 0(r), 
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In Eq. (JBJ), the derivative is taken under the constraint of constant total-particle number. 
In principle, if the free energy F is known, the distributions of particle density and order 
parameter in the superfluid state can be solved from the saddle point equations, Eq. (0EJ). 

At low temperatures and low rotation frequencies, when the Fermi energy is much bigger 
than rotation energy, Ep ^> hfl, the average kinetic energy per particle is proportional to 
n 2 / 3 (r) in the Thomas- Fermi approximation j^]. In the unitary limit, the interparticle dis- 
tance replaces the scattering length, and the average interaction energy is also proportional 
to n 2 / 3 (r). Therefore the free energy of the normal state is approximately given by 



d 3 r 



h 2 5 

CNn»(r) + V'(rWr) 

2m 



(7) 



where V'{v) = m[{uj']_ — £7 2 )r 2 + u z z 2 }/2 is the trapping potential reduced due to rotation, 
and is a dimensionless constant. From the experiment [3], the constant cn is estimated 
to be about 4.3. At low rotation frequencies, Ep 3> hfl, the constant cn does not depend on 
the rotation frequency Q. In general, the normal-state free energy Fjv should be temperature 
dependent. However, near the transition temperature which is much smaller than the Fermi 
temperature Q, the normal-state free energy F N is weakly dependent on the temperature 
and the coefficient is approximately temperature independent. 

In the unitary limit, the coefficients a and (3 in Eq. (jlj) can be expressed in terms of 
the particle density n and temperature T. The coefficient a is approximately given by 
a = c a ks{T — Tq), where c a is a dimensionless constant and ks is the Boltzmann constant. 
The coefficient (3 is a function of the fermion density, (3 = cph 2 /ftmn 1 / 3 ], where cp is a 
dimensionless constant. The transition temperature in the unitary limit is proportional to 
n 2 / 3 , ksTc = cpfi'T?! 3 /(2m), where cp is a dimensionless constant. In the experiment 0, 
the transition temperature Tc is about 25% of the Fermi temperature Tp, corresponding 
to cp = 2.4. Although the zero-temperature properties of the unitary Fermi gas has been 
numerically obtained 2^| , the accurate values of the constants c a and cp near the transition 
temperature are still unknown. In the following, we use their values in the weak-interaction 
case instead, c a = lA7cp and cp = 2.94c 2 - 24J], and discuss how these parameters affect the 
system properties near the end. In the unitary limit, the saddle-point equations (0 EJ) can 
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be further written as 

fi 2 1i p A 

- ^(V - -j— Mr) + a(r)0(r) + /3(r)|0(r)| 2 0(r) = 0, (8) 

5ft 2 ,/„ x ^ |0(r)| 2 h 2 |0(r)| 4 

bm dm rr/^r) 12m n*' 6 {r) 

where fi is the chemical potential. 

The angular momentum can be obtained from the free energy, L z = —dF/dVt. In the 
unitary limit, it is given by 

r 2<°A 
L z = J (i 3 r[mflr>(r) + 2(p*(r)(z x r ± ) • (-ifiV )<f>{r)] (10) 

= [ d 3 r{mnr 2 ± [n(r) - 4|0(r)| 2 ] - 2zft0*(r) — 0(r)}, (11) 

where 6 1 is the azimuthal angle around z-axis. When there is no superfluid current, the 
phase of the order parameter is uniform and the angular momentum totally comes from the 
normal fluid, 

L z = J d 3 rmttr 2 ± [n(r) - 4|0(r) | 2 ]. (12) 

Therefore we identify n n (r) = n(r) — 4|0(r)| 2 as the density of the normal fluid and 
n s (r) = 4|0(r)| 2 as the superfluid density In rotating systems, if the density of the an- 
gular momentum can be measured experimentally, the superfluid and normal fluid densities 
can be obtained from Eq. (jll|) . 

III. VORTEX STRUCTURE AND THERMODYNAMIC AL PROPERTIES 

In a single-vortex state, the order parameter takes the form <f)(r) = f(r±, z)e ld , where 
for convenience / can be chosen as a positive function, / = |0|. The order parameter and 
particle density of the vortex state can be obtained from the saddle-point equations (jHJ EJ) 
which can now be further simplified, 

- ^(-|^ + S + |( - 4) + 2^V1/ - 2hnf + af + Pf = 0, (13) 
2m r_i_ or± otj_ oz a rj_ 

— CN n 2 / 3 + V' - —c T c — c fA =0 (14) 

bm N ^ 3m T "n 1 / 3 12m ^n 4 / 3 

The angular momentum in the vortex state can be obtained from Eq. (fTTj) . 

L z = J d 3 r[mnrl(n - 4/ 2 ) + 2ft/ 2 ]. (15) 




FIG. 1: The density profile of the vortex at T = 0.7Tq and z = 0. The solid, dashed, dotted, and 
dotted dash lines are density profiles at rotating frequency 0Hz, 10Hz, 20Hz, and 30Hz. 

From Eq. I|15|L we can see that in a single- vortex state each supernuid particle contributes 
h/2 to the total angular momentum. 

To solve the nonlinear differential equations Eq. (JT3J) and (fT%|) . it is necessary to use 
the numerical method. Here we take the experimental parameters from Ref. ^(|, in which 
the average number of Li 6 atoms is about two million, N « 2 x 10 6 , the axial and radial 
trap frequencies are given by v z = 23Hz and v r = 57Hz. The nonlinear partial differential 
equations |T3|) and ((H}) are solved by the iteration method. 

In Fig. (JTJ), the density profiles of the vortex are plotted at T = 0.7T C and different 
rotation frequencies. The vortex produces a density dip at the center and the normal fluid 
is dominant inside the vortex core. This is because the supernuid density vanishes at the 
center, as shown in Fig. (J2J). At very low rotation frequencies, from Eq. (|13|) the size of 
the vortex core is given by £ = %/<J— 2ma(0). When the rotation frequency increases, the 
density decreases near the center of the trap as shown in Fig. (JT|). due to the increase of 
the centrifugal force. When the density decreases, the coefficients a and (3 in the Ginzburg- 
Landau theory increase, and from Eq. (|13j) the size of the order parameter decreases. As a 
result, the supernuid density also decreases with the increase of the rotation frequency, and 
the size of the vortex core increases, as shown in Fig. (J2J). 

When the temperature increases, the supernuid density decreases as shown in Fig. Q. 
The size of the vortex core becomes larger because the coefficient a increases with the 
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FIG. 2: The superfluid density profile at T = 0.7Tq and z = 0. The solid, dashed, dotted, and 
dotted dash lines are superfluid density profiles at rotating frequency 0Hz, 10Hz, 20Hz, and 30Hz. 
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FIG. 3: The superfluid density profile at Q/(2ir) = 30Hz and z = 0. The solid, dashed, and dotted 
lines are superfluid density profiles at 0.7Tc, 0.8Tc, and 0.9Tc- 

increase of temperature. The density close to the vortex core also decreases as shown in Fig. 
(J3J). To keep the total particle number constant, the densities at other places increase. As 
a result, the density at the center slightly increases with the increase of temperature. 

The total free energy as a function of the rotation frequency is plotted in Fig. (j3J). As the 
rotation frequency increases, the free energy decreases mainly due to the reduction of the 
trapping potential. As the temperature increases, the free energy also increases. In contrast, 
the absolute value of the condensation energy decreases with the increase of temperature as 
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FIG. 4: The density profile of the vortex at Q/(2ir) = 30Hz and z = 0. The solid, dashed, and 
dotted lines are the density profiles at 0.7Tc, 0.8Tc, and 0.9Tc- 
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FIG. 5: The total free energy as a function of rotation frequency. The solid, dashed, and dotted 
lines are free energies at 0.7Tq, 0.8Tc, and 0.9Tc- 

shown in Fig. @. The absolute value of the condensation energy also decreases with the 
increase of rotation frequency due to the decrease of the total superfluid fraction as shown in 
Fig. (J7J). The total angular momentum almost increases linearly with the rotation frequency 
as shown in Fig. (jHJ). The linear dependence comes from the normal fluid contribution, and 
in contrast the superfluid in the vortex state contributes a small constant to the total angular 
momentum. 
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FIG. 6: The condensation energy as a function of rotation frequency. The solid, dashed, and dotted 
lines are condensation energies at 0.7Tc, 0.8Tc, and 0.9Tc. 
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FIG. 7: The superfluid fraction as a function of rotation frequency. The solid, dashed, and dotted 
lines are the ratios of the superfluid atom number to total atom number at 0.7Tq, 0.8Tq, and 
0.9T C . 

IV. DISCUSSION AND CONCLUSION 

So far we use the weak-interaction expressions of the parameters c a and Cp to obtain vortex 
properties. It is likely that the accurate expressions of these parameters in the unitary limit 
are different. Here we examine how sensitive the vortex properties are to these parameters. 
In Fig. ©, the density profiles of the vortex at T = 0.9T C with various values of c a and cp 
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FIG. 8: The total angular momentum as a function of rotation frequency. The solid, dashed, and 
dotted lines are the angular momentum at 0.7Tc, 0.8Tq, and 0.9Tc. 
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FIG. 9: The density profile at T = 0.9Tc and £l/(2ir) = 10Hz. The solid line is the density profile 
with the weak-interaction expressions of c a and c«. The dotted-dash and double-dotted-dashed 
lines are the density profiles with c a decreased and increased by 20% respectively. The dashed and 
dotted lines are the density profiles with cp decreased and increased by 20% respectively. 

are plotted. As shown in Fig. Q, when c a increases or cp decreases, the dip at the center of 
the vortex becomes bigger. This is due to the increase of order parameter and the increase of 
the absolute value of the condensation energy as shown in Fig. (IK)jl . If the vortex structure 
can be measured in the experiment, the values of the parameters c a and c@ can be pinned 
down, which is very helpful to obtain other properties of the system. 
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FIG. 10: The condensation energy at T = 0.9Tc as a function of Vt. The solid line is the condensa- 
tion energy with the weak- interaction expressions of c a and cp. The dotted-dash and double-dotted- 
dashed lines are the condensation energies with c a decreased and increased by 20% respectively. 
The dashed and dotted lines are the condensation energies with cp decreased and increased by 20% 
respectively. 

In summary, we have developed a phenomenological theory to study the vortex properties 
of a superfluid Fermi gas in the unitary region. This charge-neutral system under rotation 
is mapped onto an electron system in a vector potential. In the framework of the Ginzburg- 
Landau theory, the order-parameter and density distributions in the equilibrium state are 
solved numerically from the saddle-point equations of the free energy. We identify that 
the superfluid density is four times the magnitude of the order parameter squared, which 
can be determined from the angular-momentum density. The size of the vortex core is 
found to increase with the rotation frequency due to the density decrease near the core. 
Thermo dynamical properties of the vortex, such as the free energy and angular momentum, 
are also obtained. Our theoretical results can be directly tested in experiments. 

Although the vortex properties in the unitary region are studied extensively in this pa- 
per, more work are needed to study the rotation effect in other regions of the BEC-BCS 
crossover, such as the crossover region and the BEC side. It is also interesting to extend 
this phenomenological approach to study the properties of vortex lattices in the BEC-BCS 
crossover. 
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